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Abstract 

In our continuing effort to understand hadronic matter at high density, we have developed 
a unified field theoretic formalism for dense skyrmion matter using a single Lagrangian to de- 
scribe simultaneously both matter and meson fluctuations and studied in-medium properties of 
hadrons. Dropping the quartic Skyrme term, we incorporate into our previous Lagrangian the 
vector mesons p and u in a form which is consistent with the symmetries of QCD. The results 
that we have obtained, reported here, expose a hitherto unsuspected puzzle associated with the 
role the u meson plays at short distance. Since the u meson couples to baryon density, it leads 
to a pseudo-gap scenario for the chiral symmetry phase transition, which is at variance with 
standard scenario of QCD at the phase transition. We find that in the presence of the uj mesons, 
the scale-anomaly dilaton field is prevented from developing a vanishing vacuum expectation 
value at the chiral restoration, as a consequence of which the in-medium pion decay constant 
does not vanish. This seems to indicate that the u degree of freedom obstructs the "vector 
manifestation" which is considered to be a generic feature of effective field theories matched to 
QCD. 

Pacs: 12.39-x, 13.60.Hb, 14.65-q, 14.70Dj 
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1 Introduction 



The recent experimental and theoretical developments in dense matter physics have shown that 
the phase diagram of hadronic matter is rich and highly non-trivial. At high temperature and/or 
density, hadrons have quite different properties than at normal conditions. Chiral symmetry is 
believed to be restored and therefore the quark condensate (qq) of QCD, its order parameter, is 
expected to drop as matter is heated and/or compressed. 

In trying to understand what happens to hadrons under extreme conditions, it is necessary 
that the theory adopted for the description be consistent with QCD. In terms of effective theories 
this means that they should match to QCD at a scale close to the chiral scale A x ~ 4:irf w ~ 1 
GeV. It has been shown that this matching can be effectuated in the framework of hidden 
local symmetry and leads to what is called "vector manifestation (VM)" p] [2] which provides a 
theoretical support to the in-medium behavior of hadrons predicted in 1991 011]. 

We have carried out several studies in addressing the problem of dense matter within the 
general scheme we are adopting. We have shown that the simple skyrmion Lagrangian can 
describe both infinite matter and pionic fluctuations thereon and their behavior as the density 
increases 0. However in this first study we found a puzzling feature, namely that the Wigner 
phase represented by half-skyrmion matter supported a non- vanishing pion decay constant. This 
was interpreted as a possible signal for a pseudo-gap scenario which according to jTj would be 
at variance with QCD and therefore could not be realized in nature. 

In [Hill], we incorporated into the standard skyrmion model the scale anomaly of QCD. In 
fact it was this feature which led to BR-scaling in 3 and therefore we were able to reformulate 
this phenomenon in a more accurate way. In this scenario /„- vanishes at the chiral phase 
transition and other properties of the vector manifestation scenario are reproduced. Thus large 
N c and scale anomaly seem to be the minimal ingredients to satisfy the matching between the 
effective theories and QCD. 

The purpose of the present paper is to do away with the ad hoc Skyrme quartic term and 
extend the model by incorporating the lowest-lying vector mesons, namely the p and the to. It 
is known that these vector mesons play a crucial role in stabilizing the single nucleon system [S] 
as well as in the saturation of normal nuclear matter [2] and moreover it would be important to 
determine how their properties behave in the medium, since they lead to measurable quantities 
that can be confronted by different descriptions. 

We consider a skyrmion-type Lagrangian with vector mesons possessing hidden local gauge 
symmetry llf)| II lj , spontaneously broken chiral symmetry and scale symmetry. Such a theory 
might be considered as a better approximation to reality than the extreme large N c approxi- 
mation to QCD represented by the Skyrme model. In our approach a skyrmion description for 
the multi-baryon system (including infinite matter) can be obtained with the parameters of the 
theory fixed by meson dynamics. Our interest here lies in the investigation of how the vector 
mesons affect the description of nuclear matter, i.e., the density and the properties of the chiral 
restoration phase transition. For this purpose we look at the ground state of the infinite matter 
system which is a soliton solution (i.e., skyrmion matter) of the Lagrangian. As one varies the 
density of the system the parameters of the theory involving this process must adapt to the 
density of the skyrmion background. We will describe their change with density in our model 
for nuclear matter. 

The content of this paper is as follows. In Section 2, we write down our model Lagrangian 
which is the simplest form of skyrmion Lagrangian, which incorporates the vector mesons in a 
way consistent with the scale anomaly and symmetries of QCD . In Section 3, a single skyrmion 
is analyzed to define the parameters of the theory at zero density. In Section 4 we analyze the 
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two skyrmion case in order to assess the behavior of the vector-meson components of the B / 1 
skyrmion. The low density realization of the skyrmion matter in our model, an FCC crystal, is 
described in Section 5. Section 6 is devoted to the study of the chiral restoration phase transition 
and the important role of the uo in the model. Some concluding remarks are given in Section 7. 



2 Model Lagrangian 

The starting point of our work is the skyrmion Lagrangian introduced in our previous work 
[5] which contains the proper, but minimal, realization of spontaneous symmetry and scale 
symmetry breaking, fundamental properties of QCD, in the effective mesonic degrees of freedom. 
To it we incorporate the vector mesons, p and u, maintaining the adequate implementation of 
the symmetry realizations. Specifically, the model Lagrangian, which we investigate, is given by 



f 2 / \ f 2 2 / \ 

1 " ' X ' Trid^d^U) + Tr(U + U ] - 2) 



4 \f x 
fl 1 



fx 



X 



4 a [j^j Tr[^ + r M + i( 5 /2)(f 4 



2 — in . 7P V — I 



( X // X ) 4 (ln(x// X ) 



/V ■ P 



4' 4 



UJ.n.UJ 



where 



U = e*v(iT -tt/U) =f, (Pa) 

^ = £ f ^£, and r„ = £8^, (JUb) 

Pnv = d^pu - d u p^ + gp^ x p u , JTJc) 

= d^u u - dvUJ^, (JUd) 

= -^e^ aP TT{U^d v UU^d a UU^dpU). JUe) 



Note that the Skyrme quartic term is not present in the model. The vector mesons, p and uj, 
are incorporated as dynamical gauge bosons for the local hidden gauge symmetry of the non- 
linear sigma model Lagrangian and the dilaton field \ ls introduced so that the Lagrangian has 
the same scaling behavior as QCD. The physical parameters appearing in the Lagrangian are 
summarized in Table. 1. Throughout this paper, we take the following convention for the indices: 
(i) a,b, ■■■ = 1,2,3 (Euclidean metric) for the isovector fields; (ii) i,j, ■■■ = 1,2,3 (Euclidean 
metric) for the spatial components of normal vectors; (iii) p,,u,--- = 0,1,2,3 (Minkowskian 
metric) for the space-time 4- vectors; (iv) a, /3, ■ ■ ■ = 0, 1, 2, 3 (Euclidean metric) for isoscalar(0)+ 
isovectors( 1,2,3). 

Let us analyze the free space meson Lagrangian, i.e. the zero baryon number (B = 0) sector. 
The vacuum values of the fields are given by, 

U vac = 1, P/i,vac = w ^,vac = 0) Xvac = fx- (~) 

Letting the fields fluctuate in this vacuum through the Ansatz : 

U <= exp(ir a vf a ), p a ^^=p a ^ u^^u^, % <= fx + X- ( 3 ) 
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Table 1: Parameters of the model Lagrangian 



notation 


physical meaning 


value 


u 


pion decay constant 


93 MeV 


fx 


X decay constant 


210 MeV 


9 


pirn coupling constant 


5.85* 




pion mass 


140 MeV 


m x 


X meass 


720 MeV 


my 


vector meson masses 


770 MeVt 


a 


vector meson dominance 


2 



* obtained by using the KSFR relation my = rtvL = = af%g 2 with 



a = 2. cf. g P TT-n = 6.11 from the decay width of p — > tttt. 
t experimentally measured values are m p —768 MeV and m^=782 MeV. 



we create the mesons as can be seen by plugging the Ansatz (JHJ) into the Lagrangian Q and 
expanding in fields and derivatives of the fields. We can then write the Lagrangian as 



where 



£ — Aree + Ant (4) 
£ free = \d^ a 3^ a - \m\^ a - \d ^ d» ~p™ + \m\p~ v a ~p a * 



and 



£ i nt = ge abc p a ^ b d^ c + ---. gb) 

This Lagrangian provides physical meaning to the parameters of the model as listed in Table 1. 
The term chosen to appear, as an example, in Eq.(j3Ja) determines the p — > 7T7t decay width as 

T-l _ ^ gpTTTT |9m| 3 /rt 

1 p— >7T7T „ a O ' V / 

3 47T nip 

where \q m \ being the momentum of the pions in the decay rest frame. 

3 The B=l Skyrmion : Hedgehog Ansatz 

The solitons of these effective theories are the skyrmions. From Eq.Q the spherically symmetric 
hedgehog Ansatz for the B = 1 soliton solution of the standard Skyrme model can be generalized 
to 

U B=1 = exp(tf -rF(r)), (6) 
p^r l =e^f^, ^ = 0, (7) 
ul=l = 0, u,^ 1 = UW(r), (8) 

and 

X B=l = f x C(r). (9) 
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The boundary conditions that the profile functions satisfy at infinity are 

F(oo) = G(oo) = W(oo) = 0, C(oo) = 1. 



(10) 



The Ansatze for p and ui can be inferred from their equations of motion by ignoring the 
space-time derivatives and taking their masses to infinity. Thus for the p, from C np , we get 



~ -Tr 



(J 



ink-k^ - cosF(r) 







(/i = 0), 



(11) 



and for the u>, from C w + Cy/Zi we have 



3/7 







^ ~ -^-£^ = { 3g sin 2 F{r) dF 
mf 



(/* = o). 



(12) 



m, 2 , 27r 2 r 2 dr 

Equations ()11I12|) . together with the behavior of F(r) near the origin, 

F(r) = ir — ar — ^r 3 + • • • , (r < 1) 
provide us with the boundary conditions for G{r) and W(r) for small r : 

G(0) = -2, W(0) = 0. 



(13) 



(14) 



The equations of motion for the various profile functions F(r), G(r), W{r) and C(r) can be 
gotten by minimizing the soliton mass : 



TP — TP B=1 _L TP B=1 _L TP B=1 _l_ f?- B=1 _L TP B=1 _l_ F B=1 



where 



B=l 



I- / dr^- \C 2 [F' 2 + A 1 sin 2 F ) + 2m 2 C 3 (l - cos f ) 

TOO 



("OO 

E^ 1 = 4n J dr2f 2 C 2 (G + 1 - cos F) 2 , 



E B = X = 4vr / dr^, { G' A + 



/2 



o 5 



G 2 (G + 2) 
2T 2 



=1 = -4vr J r 2 dr-^- \w' 2 + m 2 C 2 ^ 2 } 



E 



B=l 
WZ 



4vr / r z dr^UW 



35 , w^ 2 ^, 



= 4tt j™r 2 dr f -f {C 2 + im 2 (C 4 (lnC - i) + |)} 



(15) 

JTSJa) 
C3b) 

CUc) 

CSId) 
(Jpoje) 

C3f) 



They are, 



F" 



2 2C"\ , , 
- + — F' + m^CsinF 
r C J 

+J_ ( 4 ( G + l) si nF - sin2F) - - ^ t ^/ w'^, 



C 2- 



(16) 
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Table 2: Single skyrmion mass and various contributions to it. 



Model 


(r 2 ) 


E B=1 


Et 1 


TTp 


E B=V 


E^ 


&wz 


Ey 


7T/9-model 


0.27 


1054.6 


400.2 + 9.2 


110.4 


534.9 


0.0 


0.0 


0.0 


7r px-ioaodel 


0.19 


906.5 


103.1 + 1.4 


155.1 


504.1 


0.0 


0.0 


142.8 


7rpu;-model 


0.49 


1469.0 


767.6 + 39.9 


33.2 


370.7 


-257.6 


515.1 


0.0 


tt pu>x-^odel 


0.51 


1408.3 


646.0 + 29.2 


34.9 


355.7 


-278.3 


556.7 


64.2 



G" = m 2 p C\G + 1 — cos F) + G{G + ^ + 2) , (17) 

W" = --W' + mlC 2 W- A ll t —^J- F 1 , (18) 
2 



r " 4TT 2 f n r 2 



C" = --C + mlC 3 lnC 

rp A. 

f2 



+4 \ F ' 2 + 4 sin2 F + 3m 2 C7(l - cos F) 
fx I r 

+^(G + l-cosF) 2 -mlW 2 ^C. (19) 

Note that the contributions of uj to the mass, E^ =l and Ey^% , satisfy a virial theorem [S]. 
Equation ijTSl d) can be expressed as 

poo f 2 f 2 1 

E^= l = -4tt J r 2 dr^-W l-W" - -W + m 2 C 2 W | , 

which implies that, if W(r) satisfies the equation of motion (|18j) . 

-2E*=*=Ejj&. (20) 

The numerical results on the properties of the B = 1 hedgehog skyrmion (the mean square 
radius and the mass) are listed in Table 2 and the corresponding profile functions are given in 
Fig. 1. It is interesting to see the roles the vector mesons and the dilaton field play in describing 
the skyrmion. The uj meson provides a strong repulsion and hence makes the soliton heavier 
and bigger in size. Specifically, comparing the 7r/3 model with the itpuj model, we can see that 
when the presence of the uj increases the mass by more than 415 MeV and the size, i.e. (r 2 ), by 
more than .28 fm 2 . 

How does the dilaton affect this calculation? The irp model with much smaller skyrmion 
has a larger baryon density near the origin and this affects the dilaton, significantly changing 
its mean-field value from its vacuum one. The net effect of the dilaton mean field on the mass 
is a reduction of ~ 150 MeV, whereas for the irpoj model it is only of 50 MeV. The details can 
be seen in Table 2. The effect on the soliton size is, however, different: while the dilaton in the 
np model produces an additional localization of the baryon charge and hence reduces (r 2 ) from 
.21 fm 2 to .19 fm 2 , in the npuj model, on the contrary, the dilaton produces a delocalization 
and increases (r 2 ) from .49 fm 2 to .51 fm 2 . 

In sum, we note that when the uj is present, it plays a major role in the skyrmion properties 
with the dilaton field playing a minor role. 



5 




Figure 1: Profile functions - F(r), G(r), W(r) and C(r). 



4 The B=2 Skyrmion : Product Ansatz 

In the original Skyrme model with only pion fields, we can obtain a B = 2 configuration 
by simply taking the product of two B = 1 hedgehog U fields. This configuration can be 
a good approximation to the true solution when the two skyrmions are sufficiently far apart. 
Furthermore the energy of this configuration becomes the lowest when one of the skyrmions is 
rotated in isospin space by an angle tt with respect to the axis perpendicular to the line joining 
two skyrmions. Here we have to generalize this feature to incorporate the p, to and x fields in 
the scheme. When two skyrmions are sufficiently far apart, our Ansatz should describe each 
individual skyrmion correctly. Recalling the field configurations for a single skyrmion at infinity, 
the most natural Ansatz is 

Ub=2 <= C/s=i(n) * Ub=i(t2), 
^, B=2 ^ B=1(rl)+ B=1( _, 2)) 

-r 2 ^-rvi)+< =i (r 2 ), (21) 

where r\2 stand for the position of the centers of each skyrmion. The boundary conditions for 
B = 1 skyrmion imply that a multiplicative rule must apply to U and x> since their vacuum 
values at infinity are 1, and an additive rule to p and u, since theirs vanish at infinity. 
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0.5 1.0 1.5 2.0 
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Figure 2: Binding energy of the two-skyrmion system. Configuration A represents the case with 
no relative rotation. Configuration B represents the case where the skyrmion is rotated by an 
angle ir/2 with respect to the axis parallel to the line joining the two skyrmions. In configuration 
C, the angle is tt and the rotation axis is perpendicular to the line joining the two skyrmions. 

The relative orientation takes place in isospin space. Thus the isoscalar fields, to and x, do 
not undergo rotation. However the isovector fields tt and p could require non-trivial relative 
orientations. 

U B =2 = U B =x{ri) * AU B =i{r 2 )A\ 

(22) 

r ■ PT 2 = t ■ f=\n) + At ■ ft B =\r 2 )A\ 

where A is an SU(2) matrix. This Ansatz is a drastically simplified one and will not produce 
the lowest energy configuration of the B = 2 system, known to be of the toroidal shape |13| . 
However for the purpose of our calculation, we need only the Ansatz defined by Eqs. (|21jl and 
(|22|). which will be used to define the starting configuration of the skyrmion matter [3]. As in 
the original Skyrme model, the lowest energy for this type of Ansatz can be obtained in the 
configuration where two skyrmions are relatively rotated in isospin space by an angle tt with 
respect to an axis perpendicular to the line joining them. (See Fig. 2.) 

In Fig. 3, we show the density profile of the lowest energy configuration, in this approximation, 
which is of dumbbell shape with well-separated skyrmions. 

5 Skyrmion Matter : an FCC skyrmion crystal 

We have learned from the study on two-skyrmion system in Sec. 4 that the lowest-energy con- 
figuration is obtained when one of the skyrmions is rotated in isospin space with respect to the 
other by an angle tt about an axis perpendicular to the line joining the two. If we generalize 
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Figure 3: The lowest energy configuration obtained under the naive product Ansatz l|21|). 



Table 3: Symmetries of the FCC skyrmion crystal 





reflection 
(yz-plane) 


3-fold axis 
rotation 


4-fold axis 
(z-axis) rot. 


translation 


(x,y,z) 


* (-x,y,z) 


(y,z,x) 




(x + L,y + L,z) 


U = a + IT ■ TT — 


> (cr, — 7ri,7r2,7r 3 ) 


(cr, 7r 2 ,7r 3 ,7ri) 


((7, 7Tl, 7T 3 , -7T 2 ) 


(cr, — 7Ti, -7r 2 ,7r 3 ) 


Pi = ^aipPp ~ 


* (~PI,P2,P3) 


(P2,P3,Pl) 


(Pl,P3, -P2) 


("Pi, -P2,P3) 













this Ansatz to many-skyrmion matter, we obtain that the configuration at the classical level 
for a given baryon number density is an FCC crystal where the nearest neighbor skyrmions are 
arranged to have the attractive relative orientations 1 . We next proceed to incorporate the 
vector mesons by generalizing Kugler's Fourier series expansion method ^H]) developed for the 
original Skyrme model, to our 7T/9o;x-model. The symmetries of such FCC crystal configuration 
are summarized in Table 3. As for the /9-fields, we find it more convenient to give the symme- 
tries for the dual vector fields p p as defined in Table 3. Note: (i) p a (a = 1,2,3) has the same 
symmetries as 7r a ; (ii) the isoscalar fields cr, ujq and x share the same symmetries. One can easily 
check the symmetries for pf (or p a ) and luq by referring to the useful relations (fTTj) and (|T2|). 

An early discussion on the FCC crystal structure of nuclei can be found in I14| . See also |15| for a recent 
discussion. A quark-model description - which is similar in spirit to our approach - in which the gluonic mean 
field represented by a confining potential replaces the skyrmion field distribution is given by Goldman et al |16| 
and Kim Maltman et al |17|. We are grateful to Terry Goldman for bringing our attention to these references. 
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5.1 Fourier series expansion 



We next generalize our work in ref . ( [5J [H] ) , following Kugler and Shtrikman|18|. to the problem 
at hand by indicating the pertinent expansions for all the fields. 

i) The pion expansion: 

We obtain the pion fields 2 ^ from the un-normalized fields (ft^ which are expanded in 
Fourier series as 



4>l = 
4>l = 
4>l = 

and thereafter normalized 



flabc cos(7rax/L) cos(7t6t/ / L) cos(ircz/L), 



abc 



sin(irhx/L) cos (ir ky / L) cos(irlz/L), 



hkl 



cos(irlx / L) sm.{ixhy / L) cos(ixkz / L), 



(23) 



hkl 



a^ki cos(7rfcx / L) cos(irly / L) sin^nhz / L) , 



hkl 



£\ (, i — 



(24) 



In order to be consistent with the symmetry properties listed in Table 3, (a, b, c) should 
be all even or all odd integers, and (k,l) should be all even(odd) if h is odd(even). Fur- 
thermore, to provide a correct topological structure to the configuration, the expansion 
coefficients f3 a b c must satisfy the constraint, 

E Pabc = 0- (25) 



ii) The p expansion : 

We introduce 4> p a (a = 0,1,2,3) and using Eq. (fTT|) . we may express pf in terms of these 
fields 

pi = e abc 4> p b dd p c (l ~ = Eabc^di^Jil + 0g), (26) 



where (j> p = (p p / 'yJ2b=i(^b) 2 = fia/y^ ~ (^o) 2 - Similarly to the pion fields, (\> p a is obtained 
from the un-normalized fields § p a whose expansion in Fourier series is 



} j P P bc cos(jrax / L) cos(7t6t/ / L) cos(ttcz/L), 

abc 

E &hki siri(7r/ix / L) cos (rrky / L) cos (ttIz / L) , 

hkl (27) 
2_ . ct p hkl cos(irlx I L) sin(7r/iy / L) cos(7rkz/L), 

hkl 

2_ . c^hki cos(irkx/L) cos(irly / L) sin(7rhz/L), 



hkl 



which we then normalize. The procedure is analogous to that of the pion fields, although 
the expansion coefficients are different. 



2 We use = 0, 1, 2, 3) instead of a(— 4>q), ir a (~ 4>Z , a — 1, 2, 3). The superscript tt denotes that these fields 

are associated with the pion fields. We will introduce similar fields <j!>£ for the p fields. 
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iii) The uj and \ expansions : 

The isoscalar fields u>o and x have the same symmetry properties as the a and therefore 
their expansions are, 



X/f X 



22 Pabc cos(-7rax / L) cos(7r&y/L) cos(ircz / L) , 

abc 

c = (3* hr cos (irax jV) cos (irby/L) cos (ncz /L), 



(28) 



The expansion coefficients are determined such that the energy per skyrmion is minimized: 



where 



E/B = (E/B) n + (E/B) np + (E/B) p + (E/B)^ + {E/B) wz + (E/B) } 



(E/B) n = \ f d 3 x^c 2 Ud^l) 2 + 2mlc(l - <$)} 

JBox 2 I- > 



(E/B) 



Tvp — 4 



Box 



x did) 71 <bP x dj<b p 



(E/B) p = \ [ d 3 x-^ (diftdjtf - djcfcdt 
J Box og v 

(E/B) u = -\ I d?x\{{d lWo ) 2 + mlc 2 wl}, 

JBox 

(E/B) wz = \ f d 3 x^Uw Bo, 

J Box 2 

(E/B) x = \ j <?x f -f {(d lC ) 2 + \m 2 x [c\\nc - \) + $]} , 



(29) 

d3a) 
C3b) 

CSIc) 

(Eld) 
(HJe) 

GSIf) 



and 



^0 



12^ 2 



1 

2^2 



<Po 



^2 ^3 

s^i d^l 



2alg) 



Since there are no constraints on the expansion parameters for the u and x> a straightforward 
variational process of minimizing the energy fails. Note that as far as the energy of the system 
()29|) is concerned, ljq = comes out as the solution which minimizes the energy 3 . This is 
confirmed in Table 1 where we see that without the repulsive u, the single skyrmion mass 
becomes much lower. Therefore, the variational process always leads us to wq = 0, which 
corresponds to the energy per baryon of the irpx model and not of the present model. 

We realize that the u needs a careful treatment in the model that incorporates it self- 
consistently. We proceed therefore to treat the u in a more elaborate way. Instead of including 
the expansion coefficients (3^ bc into the minimization process, we fix them by solving the equation 
of motion for w{= wq/ f n ) at each step. The equation of motion derived from (|29j) becomes 



(-d'f +mlc 2 )w 



V 2 



2 b . 



(30) 



3 This solution cannot of course be the true solution of the model as it does not satisfy the equations of motion 
eq. lKjUl containing the inhomogeneous source term. 
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Note that c 2 and -Bo have the same symmetry structure. Explicitly, the expansions (|23j) and 
<HD lead to 

c 2 = ^ /?gbe cos(7rax/L) cos{irby/L) cos(ircz/L), (31) 

afec 

i?o = Pgfa. cos(-Kax I L) cos(Trby / L) cos(ircz/L), (32) 

abc 

2 r> 

with the expansion coefficients fi c ahc and (3 ab ° c determined from the configuration of and c. 
The equation of motion, therefore, reduces to a linear equation for the expansion coefficients 

Pabc e iven b y 

M abc,a'b'c'Pa'b'c' = Jf^^abc' ( 33 ) 

a'b'c' 



The matrix elements are 

M a bc,a>b>c> = {a 2 + b 2 + C 2 ) ( — ) <5 a f, c ,a'ft'c' + m 2 ^ 0a"b"c"fa',a",afb',b",bfc',c",c, (34) 



vr^ 2 



where 



S a ', a , if a" = 0, 

fa', a", a = S ^a",a, if a' = 0, 

^a,a'±a", if a' / and a" / 0. 

As formulated, our variational procedure is restricted to evolve in the space of parameters 
which satisfy the to equations of motion and therefore the "false" solution is not present in the 
scheme. 



6 Results 

We show in Figs. 4 and 5 the the numerical results of the energy per baryon E/B, (x) and (a) 
in various models as a function of the FCC lattice parameter L 4 . 

In the npx model, as the density of the system increases (L decreases), E/B changes little 
(see Fig. 4). It is close to the energy of a B = 1 skyrmion up to a density greater than 
po (L ~ 1.43). This result is easy to interpret. As we discussed in Sec. 3, the size of the 
skyrmion in this model is very small and therefore the skyrmions in the lattice will interact only 
at very high densities, high enough for their tails to overlap. 

In the absence of the uj, the dilaton field plays a dramatic role. A skyrmion matter undergoes 
an abrupt phase transition at high density at which the expectation value of the dilaton field 
vanishes (x) = 5 . 

In the npujx model, the situation changes dramatically as can be seen in Fig. 5. The reason is 
that the u provides not only a strong repulsion among the skyrmions, but somewhat surprisingly, 
also an intermediate range attraction. Note the different mass scales between Figs. 4 and 5. 
In both the irpcu and the irpu>x models, at high density, the interaction reduces E/B to 85% of 
the B = 1 skyrmion mass. This value should be compared with 94% in the np model. In the 

4 The baryon number density is given by pB = 1/2L 3 . Normal nuclear matter density po = 0.17/fm 3 corre- 
sponds to L ~ 1.43. 

5 In general, (y) = does not necessarily require (x 2 ) = 0. However, in our numerical results, (x) = always 
accompanies x — in the whole space. 
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Figure 5: E/B, (x) and (a) as a function of L in the models with the to. 



7rpX" m °del, E/B goes down to 74% of the B = 1 skyrmion mass, but in this case it is due to 
the dramatic behavior of the dilaton field. 

In the npuJX model the role of the dilaton field is suppressed. It provides a only a small attrac- 
tion at intermediate densities. Moreover, the phase transition towards its vanishing expectation 
value, (x) = 0, does not take place. Instead, its value grows at high density! 

In all the models analyzed thus far, (cr) goes to zero at high density. Thus, we have two 
different scenarios for the chiral transition. The models without the u, but with x-, tend to bring 
about a phase transition that is consistent with the "vector manifestation" scenario, with the 
pion decay constant f n vanishing. However, once the oj is present, the f n remains non-zero while 
(a) vanishes, a scenario that is reminiscent of the pseudo-gap realization. This is at variance 
with the VM scenario. 

Let us discuss in detail why and how the uj changes the VM-like behavior of the npx model 
at the phase transition into a pseudo-gap scenario. The term in the Lagrangian responsible for 
this behavior is the coupling of the u> to the topological current Bq, i.e. Eas.(|2lHe) and flUS g). 
To see this point, we write the equation of motion for the u> as 
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(-df+m*J)w = -^B { 



0, 



where we have replaced the space-dependent mass term m^c 2 of Ea.(|3U|) by an effective mass 
m* 2 . This equation can be exactly solved by means of the appropriate Green's function, 

Thus, w can be interpreted as a static potential generated by the source -Bo- The uj contribution 
to E/B can be expressed as 

(E/B)wz = \ I d s x[ tfV (^)\ (x f M - m * J *- *' ]) Bo(x>), (36) 
JBox J V 2 / 47r|x - x'\ 

and the term (E/B) u is, as mentioned previously, — h of this. Note that while the integral over 
x is defined in a single FCC cell, that over x' is not. Thus, unless it is screened, the periodic 
source Bq filling infinite space will produce an infinite potential w which leads to an infinite 
(E / B)w Z . The screening is done by the omega mass, m* , as can be seen from Ea. (|36jl . Thus 
the effective ui mass cannot vanish for the solution! Our numerical results reflect this fact: at 
high density the Bq-B$ interaction becomes large compared to any other contribution. In order 
to reduce it, x nas to increase, and thereby the effective screening mass m* ~ m LU (x) becomes 
larger. In this way we run into a phase transition where the expectation value of x does not 
vanish and therefore /„- does not vanish but instead increases. 



7 Concluding remarks 

This paper represents the natural continuation of our effort to understand the physics of nu- 
clear matter from the Skyrme model description. The (original) Skyrme model is an effective 
theory which represents the large N c limit of QCD. In some sense, it is the lowest order of 
a bosonized approach to this theory. Higher orders will be obtained by incorporating massive 
mesons. In particular, we know from phenomenology that the vector mesons play a crucial role 
in describing the N-N interaction |19j and infinite nuclear matter [Sj]. Thus going beyond our 
previous developments E] , the next natural step is to implement vector mesons in consistency 
with QCD |2()l I21j. The model we studied in this paper contains, in addition to the Goldstone 
pions, the scale dilaton associated with trace anomaly of QCD and the p and uj fields introduced 
according to the hidden local symmetry strategy. 

As in our previous work without vector mesons, the skyrmion matter possesses two phases: 
a low-density phase which we simply describe here by an FCC crystal 6 and a high-density 
phase which is described by a half-skyrmion CC crystal. In our previous work, the dilaton was 
crucial to realize the phase transition in a scenario consistent with the vector manifestation 
(VM) fixed point structure: The phase transition was signalled by the vanishing of both (a) and 
/*. However in the present model with the vector mesons, in particular with the uj, the mean x 
field cannot vanish and therefore, although (a) vanishes, /* does not and we seem to fall back 
to a pseudo-gap-type of picture which does not seem consistent with the VM and in that matter 
with the standard sigma-model scenario. The sole agent responsible for this puzzling feature is 
the uj meson. 

6 We refer the reader to our previous work for the discussion related to how we treat this unstable phase. 
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If the uj is not present, i.e., in the Trpx model, the behavior is the conventional one. The % 
field vanishes at the phase transition and chiral restoration is produced in the standard scenario. 
The p meson is basically a spectator at the classical level, producing little change with respect 
to our previously studied irx model [HI E] except that at high densities, once the p starts to 
overlap, the energy of nuclear matter increases due to its the repulsive effect at short distances. 
The densities have to be quite high since these skyrmions are very small. Since x vanishes at 
the phase transition, we recover the VM behavior, namely, /* = and m* = 0. 

The incorporation of the uj changes dramatically the scales. Since the uj produces a strong 
repulsion, the skyrmions become more massive and much bigger in size (Hj. Moreover the phase 
transition scenario changes dramatically from the previously described. In particular, \ does 
not vanish at the phase transition but it even increases its value. The mechanism of how this 
happens is simple and robust. It is that the uj coupling to the baryon density in our scheme leads 
in nuclear matter to a long range infinite interaction unless some sort of screening intervenes. 
However in medium the effective uj mass becomes m* = m u x- Thus if x decreases, the effective 
mass decreases and the screening decreases, thus the long-range interaction becomes stronger 
and ultimately will tend to dominate. In order to prevent this from happening, the background 
skyrmion counteracts so as to compensate for the increase in the baryon density interaction: It 
increases the value of the x field an d hence the effective uj mass, i.e., the screening. Since the 
in-medium value of the pion decay constant is locked to the mean-field value of x which does not 
vanish at the chiral transition, f w does not vanish, whereas the (a), whose vanishing is linked to 
the structure of the crystal, does to vanish at the critical point, signalling the phase transition 
as in our naive scenario 0. 

The qualitatively dramatic effect of the uj meson in the skyrmion structure of dense matter 
is at the same time disturbing and puzzling. In nuclear physics, uj has been a crucial element [H], 
so it makes one wonder what goes wrong at densities greater than that of normal nuclear matter. 
We have no clear answer at this point but we can think of three possibilities: (a) the minimal 
Lagrangian we use is inadequate in that short-distance physics is not fully accounted for. Since 
the uj degree of freedom accounts for nuclear interactions at short distance, higher derivatives 
and/or higher-dimension operators could be indispensable; (b) the Wilsonian matching of effec- 
tive field theories to QCD as discussed by Harada and Yamawaki j^j requires that the parameters 
of the effective theory be intrinsically density- and temperature-dependent. This would mean 
that the coupling constants - and not just the masses of the vector mesons - will have intrinsic 
dependence which may not be fully accounted for by the response to the background skyrmion 
matter that is taken into account in our model. This feature was noted already at nuclear matter 
density |22l |I] where the uj coupling was found to decrease along with the BR scaling mass of the 
uj in order to obtain the saturation and binding of nuclear matter; (c) in describing dense matter 
approaching chiral restoration in terms of explicit (as opposed to integrated-out) massive de- 
grees of freedom, the lowest-lying vector mesons may not be sufficient. It may be that the tower 
of vector mesons as suggested in the "open moose" model of Son and Stephanov [23] _ which is 
conjectured to be dual to QCD - may have to be implemented consistently with the symmetries 
of QCD. These points are being investigated and will be reported in future publications. 
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